Submitted to the Astrophysical Journal 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



O 
O 

(N 

00 
CO 

> 

m 

O 
O 
O 

Pl,: 

6 
& 

CO 



X 



ENRICHMENT OF THE INTERGALACTIC MEDIUM BY RADIATION PRESSURE DRIVEN 

DUST EFFLUX 

Anthony Aguirre,^ 6 Lars Hernquist|] Neal KatzJ^ Jeffrey Gardner]^ & David 

Weinberg| 

"Institute for Advanced Study, School of Natural Sciences, Princeton NJ 08540 
'Department of Astronomy, Harvard University 60 Garden Street, Cambridge, MA 02138 
c Department of Physics and Astronomy, University of Massachusetts, Amherst, MA 98105 
^Department of Astronomy, University of Washington, Seattle, WA 98195 
TJepartment of Astronomy, Ohio State University, Columbus, OH 43210 

Submitted to the Astrophysical Journal 

ABSTRACT 

The presence of metals in hot cluster gas and in Lya absorbers, as well as the mass-metallicity relation 
of observed galaxies, suggest that galaxies lose a significant fraction of their metals to the intergalactic 
medium (IGM). Theoretical studies of this process have concentrated on metal removal by dynamical 
processes or supernova-driven winds. Here, we investigate the enrichment of the IGM by the expulsion of 
dust grains from galaxies by radiation pressure. We use already completed cosmological simulations, to 
which we add dust assuming that most dust can reach the equilibrium point between radiation pressure 
and gravitational forces. We find that the expulsion of dust and its subsequent (partial) destruction in 
the IGM can plausibly account for the observed level of C and Si enrichment of the z = 3 IGM. At low- z, 
dust ejection and destruction could explain a substantial fraction of the metals in clusters, but it cannot 
account for all of the chemical species observed. Dust expelled by radiation pressure could give clusters a 
visual opacity of up to 0.2 — 0.5 mag in their central regions even after destruction by the hot intracluster 
medium; this value is interestingly close to limits and claimed observations of cluster extinction. We 
also comment on the implications of our results for the opacity of the general IGM. Finally, we suggest 
a possible 'hybrid' scenario in which winds expel gas and dust into galaxy halos but radiation pressure 
distributes the dust uniformly through the IGM. 

Subject headings: cosmology: theory — intergalactic medium — galaxies: abundances — dust: 
extinction 



1. INTRODUCTION 

Several independent sets of observations indicate that 
galaxies must lose a substantial fraction of the metals they 
produce during their lifetimes. First, metal lines in hot 
X-ray emitting gas in clusters and groups indicate that 
as much metal lies outside of galaxies in these objects as 
inside them (e.g., Mushotsky et al. 1996; Renzini 1997; 
Davis, Mulchaey & Mushotsky 1999; Buote 2000). Sec- 
ond, quasar absorption line studies imply that the inter- 
galactic medium (IGM) at z <J 3 is enriched to metallic- 
ity Z £ 1Q- 2 - 5 Z q (e.g., Songaila & Cowie 1996; Cowie 
& Songaila 1998; Ellison et al. 2000; Penton, Sticke & 
Schull 2000). Cosmological simulations indicate that this 
seems to require at least ~ 10% of galactic metals to be 
ejected (Aguirre et al. 2001a, b). Third, the strong posi- 
tive correlation between galaxies' masses and metallicities 
(e.g., Zaritsky, Kennicutt & Huchra 1994) is naturally ex- 
plained by the efficient escape of metals from low-mass 
galaxies (Dekel & Silk 1986; Lynden-Bell 1992). 

Most theoretical studies addressing this ubiquitous pres- 
ence of intergalactic metals have focused on the removal of 
metal enriched gas from galaxies; the gas may be removed 
by ram-pressure stripping, during dynamical encounters 
between galaxies, or as an outflow driven by supernovae 
and stellar winds. While dynamical removal undoubtedly 
occurs at some level (especially in rich clusters), it is not 
clear that it can account for the level of metallicity in the 



z = 3 IGM or the mass-metallicity (M-Z) relation of galax- 
ies (Aguirre et al. 2001a; but see Gnedin 1998). Metal 
ejection by galactic winds can explain the M-Z relation 
(winds escape low-mass galaxies more easily) and may ac- 
count for the observed level of IG enrichment (e.g., Cen 
& Ostriker 1999; Aguirre et al. 2001b), but it is unclear 
whether they can do this without overly disturbing the 
thermal or structural properties of the high-z IGM. 

A third metal removal mechanism, which has not previ- 
ously been treated in a cosmological context, is the ejection 
of dust grains by radiation pressure. As first pointed out 
by Pecker (1972) and Chiao & Wickramasinghe (1972), 
bright galaxies can exert a radiation pressure force on 
nearby grains that exceeds their gravitational attraction, 
forcing the grains into the galaxies' halos or beyond. Sub- 
sequent studies involving realistic model galaxies have con- 
firmed this idea, showing also that gas drag is insufficient 
to confine grains unless they start at small galactic scale- 
height (e.g., Ferrara et al. 1990; Shustov & Vibe 1995; 
Davies et al. 1998; Simonsen & Hannestad 1999). 

All of these studies support the idea that much of a 
galaxy's dust may be ejected during its lifetime, so it is in- 
teresting to assess the possible IG enrichment that would 
ensue. Unlike winds, enrichment by dust (partially de- 
stroyed in transit or by the IGM) would not impact the 
thermal/structural properties of the IGM or galaxies. In 
this Letter, we assess the amount and distribution of met- 
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als transferred to the IGM as dust driven by radiation pres- 
sure, using two smoothed-particle hydrodynamics (SPH) 
simulations. The first has 128 3 dark matter particles and 
128 3 SPH particles in a (17Mpc) 3 box, and ends at z = 3. 
The second, ending at z = 0, has 2 x 144 3 particles in 
a (77Mpc) 3 box. Both assume fi A = 0.6, n b = 0.047, 
Q m = 0.4, h = 0.65 and as = 0.8. The simulations are 
described in more detail in Aguirre et al. (2001a) and in 
Weinberg et al. (1999). Section || describes the method of 
adding metals and dust to the already completed simula- 
tions. Section |^ gives results pertaining to the enrichment 
of the z = 3 IGM and the z = intracluster medium, 
in several representative models. We discuss these results 
and their implications in § ^. 

2. METHOD 

The method by which we calculate IGM enrichment is 
discussed in detail in Aguirre et al. (2001a). Briefly, our 
method post-processes a limited number of outputs from 
already completed SPH cosmological simulations that in- 
clude star formation. We assume that each unit of forming 
stellar mass instantaneously generates units of metal 
mass. We then deposit this metal mass in gas particles 
near the forming star particle as follows: 

1. A fraction (1 — Y e j) of the metal is distributed in 
the nearest 32 gas particles, using the SPH smooth- 
ing kernel (see Hernquist & Katz 1989). Half of 
the locally-distributed metal is added in the form 
of dust, the other half as gaseous metal. 

2. The remaining mass is tallied for a given galaxy, 1 for 
which we also compute the mean metallicity {Z) ga \ 
and the UV-optical-NIR luminosity, using the mod- 
els of Bruzual and Chariot 2 and a Scalo or Salpeter 
initial mass function (IMF). 

3. Using (iT)gai we apply a dust correction to the lu- 
minosity from Heckman et al. (1998; see Aguirre et 
al. 2001a), normalized to give the observed ratio at 
z = in the cosmic UV-optical-NIR and FIR back- 
grounds (which are also output by the simulations). 

4. We assume a grain size distribution (GSD) in mass 
dm(a)/da and opacity (per unit mass) law from Kim, 
Martin & Hendry (1994) and Laor & Draine (1993), 
respectively, for either graphite or silicate grains. 

5. The fraction Y e j of metal formed in a galaxy is dis- 
tributed as dust spherically about the center of star 
formation. A dust mass proportional to dm(a)/da is 
placed in a shell where the radiation pressure on a 
grain of radius a balances the galaxy's gravitation. 

The process is repeated for each galaxy at each time 
step. New stars are formed with the metallicity (including 
dust) of the gas from which they form. Each gas particle 
has an accumulated mass of gaseous metals and dust, and 
we track the GSD for each particle using a 9-point piece- 
wise power law fit (see Aguirre et al. 2001a for details). 
The GSD is modified as the dust is converted to metals by 

1 B y 'g-a.la.vy' wp mpan a group of bound particles found 
hpcc.as tro. Washington. odu/took . 

2 The models are available via anonymous FTP from ftp.noao.edu. 
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Fig. 1. — Enrichment of the IGM plotted in four ways. Panel A: 
Random subsample (1 in 500) of particle metallicities for the fidu- 
cial model with graphite grains, versus overdensity <5 = p ga s/(Pgas). 
Top axis (here and in all panels) gives approximate log N(H I), us- 
ing the relation of Dave et al. (1999). The solid line shows the 
median metallicity versus S. B: As for panel A, for silicate grains. 
C: Median metallicities versus 8 for models with graphite and sili- 
cate grains, but for total (dust+gas) metal content, and for gaseous 
metals only. The shaded box roughly indicates the metallicity of 
low-column density Lye? absorbers (Ellison et al. 2000). D: As for 
panel C, but mean metallicities are plotted. E: As for panel C, but 
gives mean metallicities times the fraction of baryons at a given 8, 
showing the contribution by components with different 8 to the cos- 
mic metal density. The thick line shows the distribution assuming 
constant metallicity (with the same total metal mass). 



thermal sputtering by the IGM (using the yields of Jones 
et al. 1994), or as new (unsputtered) dust is added to the 
particle. 



3. RESULTS 

Our basic model assumes graphite grains, a 1:1 ratio be- 
tween the cosmic UV-optical-NIR and FIR backgrounds at 
z = (c.f. Madau & Pozzetti 2000), a Scalo IMF with cut- 
offs at 0.1 M Q and 100 M Q , y, = Z Q , and Y cj = 0.5. The 
last assumption is maximal, as only ~ 1/2 of a typical 

using the SKID package, publicly available at ^lttp: / /www- 



Fig. 2. — Left: Maximal dust ejection radius h^ ust vs. galaxy 
mass for 2 = 3. Right: /Jd ust vs. mean metallicity. 

galaxy's metals are in dust. We also give corresponding 
results for silicate grains. 

Figure |l| shows the key results at z = 3, using the 128 3 
simulation. The top two panels give a sparse sampling 
of individual particle metallicities, versus the gas overden- 
sity 8. The stellar yield y* is uncertain by perhaps a fac- 
tor of two, and all of the curves could be scaled vertically 
for a higher assumed value. The metallicity at 5 <J 10 4 
could also be (roughly) scaled by Y Ci for lower assumed 
values. The bar at the bottom of each panel shows the 
zero metallicity particles and indicates that the distribu- 
tion is rather inhomogeneous, especially for silicate grains 
(Panel B). This can also be seen by comparing panels C 
and D, which show the median and mean metallicity vs. 
8. The latter shows that dust ejection can provide enough 
metals to account for the Lya observations (indicated by 
the hatched rectangle), though the enrichment may not, in 
these models, be uniform enough. It is important to note, 
however, that (assuming grains decouple from the galac- 
tic gas) our method always underestimates the radius to 
which the grains can escape, because they would inevitably 
reach the force balance radius with some velocity and over- 
shoot it. Thus the distribution should probably be more 
uniform than shown here. (The ejection radius should also 
be limited by the average dust velocity Vd times the avail- 
able time, but introducing this limit does not change the 
fiducial model results unless Vd <J lOOkms -1 , slow com- 
pared to velocities seen in more detailed studies of dust 
ejection.) 

Figure [j] gives results for both the total metal enrichment 
(dashed and solid lines), and for the gas-phase enrichment 
(single- and triple-dot-dashed lines), where grains have 
been converted to gas by thermal sputtering only. Because 
destruction by both thermal and nonthermal sputtering 
during grain ejection would destroy more dust, true gas- 
phase abundances should lie above the latter two curves 
(although if grains are destroyed very efficiently at small 
radii they will not survive to pollute the low-density re- 
gions). 

The models with different dust corrections (e.g. chang- 
ing the 1:1 ratio in cosmic backgrounds to 1:2 or 3:1) give 
differences in z = 3 enrichment comparable to the 



Fig. 3. — Dust extinction for a simulation cluster at z = 
in the fiducial model, with graphite grains. Both images are 
1.2Mpcxl.2Mpc, projected through a 1.2 Mpc cube. The left sur- 
face gives visual extinction, assuming Ky = 4 X 10 4 cm 2 g — 1 . The 
right panel simulates what a sheet of white paper would look like 
through the dust of the cluster. 

differences between the silicate and graphite models 
(which differ in dust opacity by a factor of a few). Sim- 
ilar changes are induced by different assumed IMFs (see 
Aguirre et al. 2001a). 

Quantities pertainingto the galaxies ejecting dust at 
z = 3 are shown in Fig. |2[ The left panel, giving the max- 
imal dust ejection radius vs. the galaxy (baryon) mass, 
shows that ejection is most effective from the larger galax- 
ies. This, and a correlation between mass and metallic- 
ity (resulting from more efficient star formation in larger 
galaxies), largely washes out the anti-correlation between 
ejection radius and metallicity one would expect from the 
metallicity-dependent dust correction (as shown in the 
right panel). At low redshift, the balance reverses, and 
high mass galaxies eject metals slightly less efficiently due 
to their large dust corrections. 

The 144 3 simulation (which runs to z — 0) allows us 
to assess the enrichment of the low-z IGM by dust ejec- 
tion. This simulation only resolves galaxies of baryon mass 
j> 1O 1O - 7 M , but these galaxies dominate the observed 
z = mass function, and the more efficient ejection of 
grains from large galaxies at high-z (when small galaxies 
contribute relatively more mass), so we capture the bulk 
of the metal enrichment. In the fiducial model described 
above, the ICM of the most massive groups/clusters is en- 
riched to w 1/5 Zq. 

Dust is destroyed efficiently in the hot ICM, but enough 
remains that some extinction can occur. Figure |j| shows 
the optical depth through a simulation group/cluster of 
baryonic mass 2 x 10 13 M© in the fiducial graphite model, 
assuming a dust visual opacity of ny = 4 x 10 4 cm 2 g _1 
(reasonable for a more realistic silicate- graphite mixture). 
Except along paths through galaxies, the cluster opti- 
cal depth is typically < 0.2 mag; the central region has 
Ay ~ 0.5 mag. Poorer groups show less opacity. We note 
also that the same model predicts a general 'diffuse' extinc- 
tion to z = 0.5 of w 0.05 - 0.1(kv/4 x 10 4 cm 2 g" 1 ) mag, 
which is comparable to the difference between Hubble dia- 
grams for different cosmological models at z = 0.5 (Aguirre 
1999) and could potentially be important in observational 
cosmology. 3 

4. DISCUSSION AND IMPLICATIONS 



3 The far-infrared emission from such dust would not violate constraints from the observed far-infrared or microwave backgrounds; see 
Aguirre & Haiman (1999). 
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In § [l] we argued that dust ejection is an interesting 
alternative to dynamical or wind enrichment of the IGM 
because it may be efficient, yet not disturb the IGM or 
galaxies in a way incompatible with observations. Our 
simulations, which produce reasonable predictions for the 
masses, luminosities, and spatial distribution of galaxies, 
support this possibility by indicating that most galaxies 
at high z have properties that would tend to repel dust 
grains, out to a radius large enough that the low-density 
IGM can be significantly polluted. The chief uncertainties 
in our calculation are not the detailed choices of dust opac- 
ity, IMF, dust correction, cosmological parameters, etc. 
(all of which are probably uncertain only at a level which 
does not significantly affect our results), but rather in the 
physics of dust ejection itself. We assume that most dust 
reaches the point of equilibrium between radiation pres- 
sure and graviational forces, but realistically dust might 
be destroyed in transit, or confined to galaxies by other 
forces. Gas drag can confine grains at small galactic scale- 
heights, but this still allows a large dust outflow when the 
circulation of gas in the galaxy is considered (Shustov & 
Vibe 1995). But magnetic fields (not included in our treat- 
ment) might be extremely important, since charged grains 
in a microgauss field would oscillate about field lines with a 
Larmor radius significantly smaller than the galaxy scale. 
To escape, dust may diffuse along a vertical component 
of the magnetic field (Shustov & Vibe 1995), perhaps en- 
hanced by low-level winds or by Parker instabilities (Chiao 
& Wickramasinghe 1972; Ferrara et al. 1991). Magnetic 
fields may also be much weaker at high-z if they have been 
amplified by a dynamo since then. 

If dust can escape magnetic fields, our calculations show 
that it could significantly pollute the IGM. A unique sig- 
nature of enrichment by dust is that while dynamics or 
winds would pollute the IGM with chemical abundances 
similar to those of the galaxies, dust ejection can only en- 
rich the IGM with elements such as C, Si, and Mg, which 
solidify as grains. Elements such as N, Zn, and the noble 
gases, which are very lightly depleted onto grains, should 
only be ejected in trace amounts. Thus by measuring the 
relative ratio of N to C or Si in Lya lines, one could con- 
strain the pollution by dust. Presently N is detected only 
in absorbers of fairly high (N(H I) ~ 10 16 cm~ 2 ) column 
density (Songaila & Cowie 1996), but pushing these obser- 
vations to lower H I columns could give strong constraints 
on (or evidence for) dust enrichment. 4 At low redshifts the 
significant abundances of Ne and Ar in cluster gas (e.g., 
Mushotsky et al. 1996) indicates that dust cannot be the 
sole pollutant of the ICM and that some enrichment by 
other mechanisms must occur. Higher quality data from 
Chandra should allow a much more interesting test of the 
importance of dust ejection. 

Our calculations also give a fairly accurate assessment 
of the expected opacity of rich clusters if about half of the 
observed enrichment were due to dust ejection (unless the 
grain opacities are significantly higher than we have as- 
sumed). Our estimate of ~ 0.2 — 0.5 mag in the central 
few hundred kpc of clusters is roughly comparable to both 



claimed detections of cluster dust using extinction of back- 
ground quasars (e.g., Boyle et al. 1988; Romani & Maoz 
1992) or IR emission (Stickel et al. 1998), and to upper 
limits based on reddening (e.g. Maoz 1995). 5 This indi- 
cates that the general picture of substantial dust ejection 
from galaxies might provide an interesting level of extinc- 
tion through the IGM, but would not violate any current 
constraints on cluster dust density. 

The primary difficulty with dust ejection as an explana- 
tion for the Si and C enrichment of the low-density IGM 
is that it is not at all clear - theoretically or observa- 
tionally - that dust really can decouple from galactic gas; 
but if it can, there appears to be no reason why it would 
not escape to large radii. Galactic winds, on the other 
hand, are clearly observed both locally and at high-z, and 
they should certainly be able to pollute (at least) the ha- 
los of their progenitor galaxies. But spreading the metals 
to large distances may disrupt the IGM more than obser- 
vations allow (e.g., Theuns, Mo & Schaye 2000; Aguirre 
et al. 2001b). This suggests a possible 'hybrid' scenario 
in which gas and dust are expelled into a diffuse mixture 
in the halos of galaxies. But while gas remains there, the 
dust could continue, driven by radiation pressure, to large 
distances. 

For example, imagine a representative z — 5 galaxy of 
baryonic mass 5 x 10 9 M Q and UV-optical-NIR luminosity 
2.5 x 10 10 Lq driving a wind of velocity v — 200 km s _1 at 
small radii. Such a galaxy could reasonably drive a wind to 
~ 1 — 100 kpc (using the results discussed in Aguirre et al. 
2001b), but our calculations show that radiation pressure 
could exceed gravitational attraction out to 100-200 kpc. If 
0.1 /jm graphite grains were to decouple from the gas near 
the disk (say at 10 kpc), they could reach (100 — 200) kpc 
after w (0.2 - 0.4)Gyr, with velocity w 470 - 500 km s" 1 
(in calculating this we assume here that the enclosed mass 
is proportional to the radius). After ~ 1 Gyr (i.e. at 
z = 3) the grains could reach up to ~ 420 kpc; silicate 
grains could reach up to ~ 280 kpc during the same time. 
This is an upper limit since we have neglected gas drag; 
adding gas drag appropriate for 5 = 100(10) gas reduces 
the graphite distance to 250(385) kpc and the silicate dis- 
tance to 220(270) kpc. 6 but indicates that radiation pres- 
sure can quite plausibly eject dust far enough to pollute 
the IGM quite uniformly while disturbing the IGM only 
near the galaxy. 

In summary, our calculations indicate that galaxies at 
high redshift tend to repel rather than attract dust grains. 
If a substantial fraction of dust can reach at least the equi- 
librium radius between gravitational and radiation pres- 
sure forces, then the ensuing enrichment can account for 
the mean level of C and Si observed in the IGM at z ~ 3. 
Dust ejection would also enrich groups and clusters sub- 
stantially, though radiation pressure cannot account for 
all of the metals observed. The resulting dust extinction 
would be < 0.5 mag through the cores of rich clusters. 
Dust ejection and the ejection of metals by winds are, in 
some sense, complementary. Winds almost certainly drive 



4 Unfortunately (for this application), N might also be lacking if it is underproduced in the massive (perhaps low-metallicity) stars responsible 
for the enrichment at high-z; see Arnett (1995). 

5 The conclusions of reddening studies are vulnerable to changes in the dust grain-size distribution by dust destruction; see Aguirre (1999). 

6 For this estimate we use a drag force apv 2 where a and v are the dust geometrical cross section and velocity, and p is the density of the 
medium. This neglects Coulomb drag and is good in the limit of highly supersonic grains. 
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gas into the halos of galaxies, but may overly-disturb the 
IGM if the gas travels to very large radii. Dust may be 
confined to galaxies by magnetic fields or gas drag but 
should leave unimpeded if first moved into the halo. Ra- 
diation pressure acting on dust can therefore help enrich 
the IGM more uniformly than winds alone. Because only 
certain elements form dust, the possibility of intergalactic 
enrichment by dust can be robustly tested - in principle - 
by measuring ratios between refractory and non-refractory 
elements in the IGM. 
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